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Hofava-Lifshitz gravity is a potentially UV complete theory with important implications for the very early 
universe. In particular, in the presence of spatial curvature it is possible to obtain a non-singular bouncing 
cosmology. The bounce is realized as a consequence of higher order spatial curvature terms in the gravitational 
action. Here, we extend the study of linear cosmological perturbations in Hofava-Lifshitz gravity coupled to 
matter in the case when spatial curvature is present. As in the case without spatial curvature, we find that there 
is no extra dynamical degree of freedom for scalar metric perturbations. We study the evolution of fluctuations 
through the bounce and show that the solutions remain non-singular throughout. If we start with quantum 
vacuum fluctuations on sub-Hubble scales in the contracting phase, and if the contracting phase is dominated by 
pressure-less matter, then for A = 1 and in the infrared limit the perturbations at late times are scale invariant. 
Thus, Hofava-Lifshitz gravity can provide a realization of the "matter bounce" scenario of structure formation. 
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I. INTRODUCTION 

Hofava has proposed a simple model of quantum gravity I H 2] which is conservative in the sense that it is based on using the 
usual metric degrees of freedom in four space-time dimensions, but is radical in the sense that it abandons general covariance and 
local Lorentz invariance. Instead, the theory is based on a scaling symmetry in which space and time scale differently. Spatial 
diffeomorphism and space-independent time reparametrizations remain as symmetries of the theory. Hofava-Lifshitz gravity, as 
this theory is now called, has a free-field ultraviolet (UV) fixed point. It is argued that there is also an infrared fixed point in which 
the action reduces to that of General Relativity and in which local Lorentz symmetry and space-time diffeomorphism invariance 
emerge. There have been several general studies of Hofava-Lifshitz gravity Q] and a number of studies of cosmological aspects 
of the theory 0]. 

Since Hofava-Lifshitz gravity has the same number of fields as General Relativity but has a reduced symmetry, we should 
expect an extra physical mode jl]. This mode could be ghost-like JH, it could be strongly coupled or it could be simply 

well-behaved but phenomenologically ruled out. In a previous paper |@] (see also J^l]) we showed that, in the absence of spatial 
curvature, the extra mode is not propagating at all (this conclusion was later confirmed in ifioll ). Our analysis showed that the 
strong coupling instability discussed in (6] and the ghost-like evolution studied in j5t] are regulated by taking into account the 
expansion of space which is inevitable in the presence of background matter 1 . 

Earlier, it had been shown ITLzTl (see also foi l that in the presence of spatial curvature it is possible to obtain a non-singular 
bouncing cosmology. At the bounce point the expansion rate of the universe vanishes and hence the question arises as to whether 
linear cosmological fluctuations are well-behaved in a bouncing Hofava-Lifshitz cosmology in the same way as they are well- 
behaved in a spatially flat expanding cosmology. In this paper we will firstly show that the presence of spatial curvature does 
not change the conclusion that there are no extra propagating degrees of freedom. Secondly, we show that the fluctuations pass 
through the bounce smoothly in spite of the fact that a term in the equations of motion blows up. 

As argued in fl2tl , Hofava-Lifshitz gravity in the presence of non-vanishing spatial curvature may yield a concrete realization 
of the "matter bounce" (see ITil fl5l flal for original works, lfl7ll for more recent studies and lfl8tl for a short review) scenario. 
In this scenario, fluctuations which originate as quantum vacuum perturbations of a matter scalar field on sub-Hubble scales 
in a matter-dominated contracting phase will evolve into a scale-invariant spectrum of curvature perturbations at later times in 
the expanding phase, with a special shape and distinguished amplitude of the three point function lfl9ll . However, in Il2ll the 
evolution of fluctuations was considered in the context of the Einstein gravity equations, and without analyzing their propagation 
through the actual bounce. The results of the present paper show that the equations of General Relativity indeed provide an 
excellent approximation to the actual evolution for IR modes of interest to current cosmological observations. 

The outline of this paper is as follows. We first briefly review Hofava-Lifshitz gravity. Next, we analyze the conditions which 
must be satisfied in order to obtain a non-singular bouncing cosmology. We find that in order to realize a non-singular bounce, 
non-trivial spatial curvature (either a closed or an open universe) is needed. We also specify the conditions on the matter content 
in the contracting phase which must be satisfied in order to obtain a bounce. In the next section we then extend the theory of 
linear adiabatic cosmological perturbations |8[] to the case in which spatial curvature is present. We show that there is no extra 
propagating degree of freedom, as in the case studied in ]8[] 1 . However, at the bounce point some of the coefficients in the 
equations of motion blow up. Thus, in section [VJ we study the evolution of cosmological fluctuations through the bounce. We 
find that on IR scales relevant for current cosmological observations the evolution of fluctuations in the pre-bounce contracting 
phase is indistinguishable from what happens in General Relativity. Then, we show that the fluctuations evolve smoothly through 
the bounce. Finally, we show that a scale-invariant spectrum of curvature perturbations emerges in the case of a "matter bounce", 
i.e. a bouncing cosmology in which the contracting phase is dominated by pressureless matter. There are corrections of order 
A — 1, where A = 1 is the IR fixed point at which the IR part of the action reduces to that of General Relativity. Finally, we 
discuss our results and give some conclusions. 

II. BRIEF REVIEW OF HORAVA-LIFSHITZ THEORY 

In Hofava-Lifshitz gravity space and time are treated differently. The space-time manifold has an extra structure, namely 
a given foliation of space-time into constant time hypersurfaces. Instead of full space-time diffeomorphism invariance, the 
symmetry of the Hofava-Lifshitz theory is foliation-preserving diffeomorphisms, which consists of (time-dependent) spatial dif- 



1 After the work reported in this paper was completed, a paper appeared fTHl showing how the strong coupling instability can be resolved by adding extra terms 
to the original Hofava-Lifshitz action. 

2 Note that we are considering the "non-projectable" version of Hofava-Lifshitz gravity. 



3 



feomorphisms and space-independent time reparametrizations. A key ingredient in the theory is the anisotropic scaling symmetry 

t->l*t, x l ^lx l . (2.1) 

In order to obtain a power-counting renormalizable theory of gravity in four space-time dimensions we set the scaling dimension 
z = 3. In this case, the theory in the UV region should flow to a free-field fixed point and is renormalizable by power counting. 
Meanwhile, in the IR region the theory is expected to flow to the General Relativity limit where A = 1. 
The basic variables are the spatial metric g- t j, the shift vector N l and the lapse function N. 



ds 2 



(2.2) 



The spatial metric and the shift vector are functions of space and time. For the lapse function there are two choices: either N 
depends only on time (when the so-called "projectability condition" is satisfied), or it is taken to depend on both space and time 
(the general case). We will consider the general case 3 . 

The action of Horava-Lifshitz gravity contains a "kinetic" part and a "potential" part, 



- Q9 
K 



s 9 = s 



(2.3) 



The action contains the terms consistent with the symmetries of the theory (in particular spatial diffeomorphism invariance) and 
with the correct scaling dimension. The kinetic part is given by 



■ J dtd 3 x^/g~N (K %J K 13 - \K- 



(2.4) 



where 



is the extrinsic curvature and K = g y i£y. In General Relativity, general covariance requires A 
dynamical and thus runs as the energy scale changes. 

We will take the potential part of the action to be of the "detailed-balance" form 



(2.5) 

1. The coupling constant A is 



Sy = 



dtd 3 Xy/gN 



K 2 fl 2 



where CV, is the Cotton tensor defined by 



(1-3A) 



1 -4A 



-R z 



AR - 3A' 



(2.6) 



Akl 



V9 



R\ 



\ Rs i 

4 1 



(2.7) 



Here and in the above, tensors like R are understood to be constructed from the spatial metric, and g is the determinant of 
the spatial metric. The "detailed balance" condition reduces the number of terms in the potential. The most general potential 
is discussed in llzfll . Choosing the simple form of the potential will simplify our equations (which are already complicated 
enough) but will not change our basic conclusions concerning the number of dynamical degrees of freedom and concerning the 
non-singular behavior of the solutions through the bounce. 

We consider the simplest form of matter to be coupled to gravity, namely a scalar matter field 4 . The general structure of 
the action of scalar-field matter in Horava-Lifshitz theory contains two parts: a quadratic kinetic term invariant under foliation- 
preserving diffeomorphisms and a potential term: 



= / dtd 3 x ^/g~N 



1 



2N 2 



{(p - N^iip) 2 + F{tp, di<p, gij) 



(2.8) 



where F will contain higher order terms in spatial derivatives consistent with the symmetries and with power-counting renormal- 
izability. 



3 As discussed in [ 2, 20] there might be problems in the general case when attempting to quantize the theory. 

4 As is well known and as is reviewed at the beginning of Section V, a scalar field oscillating about the minimum of its potential yields a matter-dominated 
equation of state provided that the quadratic term in the expansion of the potential about the minimum does not vanish. 
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The speed of light in Hofava-Lifshitz theory can be obtained by comparing the action with that of General Relativity. The 
Einstein-Hilbert action in 3 + 1 dimnsions is written in ADM form as 

3 C 1 A 

Seh = t^t^ / cdtd 3 x^N{-(K l3 KV - A' 2 ) +R- 2^f} (2.9) 
107mG J c c z 

The expressions for the gravitational constant and the speed of light in Hofava-Lifshitz gravity can be derived by comparing the 
coefficients in the action with those in General Relativity. In the infrared limit one obtains 



which can be seen from the ratio of the coefficients of the kinetic term and the R term. In addition, 



1 - 3A 

and 



lQ wG= ^LJ~ } (2.11) 



3k V A 2 3 2A m1 „ 



32(1 -3A) 2 
Finally, it is easy to get the coefficent of the A 2 term: 

III. MATTER BOUNCE BACKGROUND 

In this section we analyze the background cosmology of Hofava-Lifshitz gravity and study under which conditions a non- 
singular bounce will occur. 

We take the background metric to be 

ds 2 = -dt 2 +g ij dx l dx J , (3.1) 

with 

Qij — & hij ~ ~ ' (3.2) 

where r 2 = 8ijX t x 3 and k is the spatial curvature which takes the values k = —1,0, 1. As we will see, in order to obtain a 
matter bounce in Hofava-Lifshitz gravity, k ^ is needed. Note that we are using units in which the spatial coordinates Xi are 
dimensionless (with respect to the usual dimensions - not the anisotropic scaling dimension) but the scale factor carries dimension 
of length. 

The background equations of motion take the form 

6(3A-1) 2 _ 3«y fk V 

k 2 9 8(3A-1) I, a 2 J ' (3.3) 

p + 3H(1 + w)p = , 

where p and p are the energy and pressure densities, respectively, and the equation of state parameter w is w = p/p. All other 
background equations are consistent with the above equations, for example 

2(3A-1) ■ _ (l + w)p gV fk \k 

~^>— H — + 4(3A^T) U " ) ~ 2 ' ( } 

From (13.31 ) it follows that a bounce can only occur if k ^ (since otherwise H = cannot be obtained). 
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In this work, we consider scalar field matter. The background energy p and pressure p for this matter take the form 

•2 -2 

P=^- + V, p=^--V. (3.5) 

Since the cosmological constant must be tuned to be very small today, we will concentrate on the case when k/a 2 A is always 
satisfied. 

When the equation of state for the scalar field satisfies w < 1/3, then in the contracting phase the higher order curvature term 
in ( 13.3b will eventually catch up with the energy density p, resulting in a cosmological bounce, a time when H = and H > 0. 

To take one step further, we would not like to have super-deflation or super-inflation around the bounce. For this purpose, we 
need H to change sign twice, once before and once after the bounce time when H = 0. This is achieved if 

k ,\ 2 4 ( k . \ k 



A < ^n ^ ~~ A - (3 - 6) 

i J 3(1 + w) \a z J a z 

which for negligible cosmological constant is realized is w < 1/3. Otherwise, the cosmology will either begin with a phase of 
super deflation leading to a bounce and then to deceleration, or with accelerated contraction followed by a bounce and then super 
inflation. Note that the "matter bounce" conditions k/a 2 ^> A and w = yield a usual bounce without super deflation or super 
inflation. 

There exist three different phases in a matter bounce cosmology: the first is the contracting phase during which the equation 
of state is dominated by pressure-less matter. This phase ends when the spatial curvature-induced higher derivative terms in the 
equations of motion become important. When this occurs, the second phase - the bouncing phase - begins during which the 
curvature-induced terms will allow the universe to evolve from contraction to expansion in a non-singular way. The last phase 
begins when the curvature-induced higher derivative terms cease to be important as the universe grows in size. At that point, the 
expanding phase that we observe today begins. In order to be consistent with late time cosmology, there needs to be entropy 
generation during or after the bounce such that we get an expanding radiation phase. How to generate the required amount of 
entropy is an issue we will not address here. 

In the contracting and expanding phases, the scale factor can be parameterized as a power law: 

a = a B i]^ , (3.7) 

which yields 

^ (1 + 3w)(r] - rf B ) ' < * 3 ' 8 ' ) 

where rfs is the time of the bounce. Note that for a matter-dominated contracting phase the equation of state parameter is w = 0. 

As in |[l7ll . we model the evolution of the Hubble parameter in the bouncing phase by linearly expanding in time about the 
bounce point: 

= 1 ~v ( 3 - 9 > 

l-y{i]-ri B y 

which leads to 

H= * y(r ;- riB \ 2 (3.10) 
l-y[ri - f]B) 

In the following sections we will study the evolution of linear cosmological fluctuations in this background. We will assume 
that the bounce occurs at a radius which is large in Planck units (this is a natural assumption if the universe starts out cold and 
with a length scale related to the initial temperature by dimensional analysis). Later on in the text we will call this the "large 
bounce radius assumption". 

IV. PERTURBATIONS IN THE PRESENCE OF CURVATURE 

We will focus on scalar metric perturbations. In General Relativity, these fluctuations can be described in terms of four scalar 
functions of space and time <fi, tp, B and E (see e.g. l22tl for a comprehensive review of the theory of cosmological perturbations 
and i23tl for a shorter overview): 

ds 2 = -(1 + 2cj})dt 2 + 2V i Ba(t) 2 dtdx i + a{tf [(1 + 2^)5 l3 + 2V l V J E]dx i dx J . (4.1) 
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There are two scalar gauge degrees of freedom which allow us to elimate two of these four functions. For example, in longitudinal 
gauge one chooses to set E = B = 0. However, in Hof ava-Lifshitz gravity one loses one of the gauge degrees of freedom, namely 
the one corresponding to space-dependent time reparametrizations. Thus, one can only eliminate one of the scalar degrees of 
freedom and one should expect an extra propagating mode. 

We will follow [8] and use the remaining gauge freedom in the scalar sector to eliminate the function E. Thus, we write the 
perturbed spatial metric in the form 

5 . y ^(l-2^=a 2 ^M^. (4.2) 
(l + £r 2 ) 

Due to the conformal properties of the Cotton tensor, for the perturbed metric ( 14.2b CV, = and e y RaV jR\ = 0. 

In addition to the fluctuation in the spatial metric, there are perturbations of the shift vector, the lapse function, and the matter 
scalar field: 

N = l + (/){t,x), 

N t = V<B(i, x) , (4.3) 
<P = <A) + Q{t,x) . 

Note that we are not enforcing the "projectability condition". If we had enforced this condition, then <fi would be constrained 
to be a function of time only, and we could use the residual gauge freedom of space-independent time reparametrizations to set 

<p = 0. 



A. Solving the Constraints 



The equations of motion for N and iVj are: 




4 (KhK* XK 2 ) ±-C, ,(-> + P^e^RaV.Ri K> 



2w 4 



2w 2 ' 



2 ,,2 



RijR 13 



;i-3A) 



1 -4A 



R + AR — 3A 2 - 



2N 2 



(<p-N%<p) 2 + F, 



For the background metric d4~2t . dj = and e ljk R a V J R l k = 0. 
At first-order, the energy constraint gives 



0= 4(l-3A)g AB . 



12ff 2 (l-3A) , . 2 
12 + ( A) 



k 4 (k - a 2 A) fi 2 (a 2 A?/' + 3fcV) - 2a 4 (-l + 3A) (l2H(-l + 3A)^ + n 2 (d<p + QV'^j 



2a 4 K 2 (-l + 3A) 



while the momentum constraint yields 

o = i 



(-1 + 3A) [H(f> + i> 



2k 



B + (1 - X)ABj 



Q ■ 



(4.4) 



(4.5) 



(4.6) 



In the above A is the Laplacian constructed using the background spatial metric gij. 

As was done in the spatially flat model in iH, we can combine the above two constraint equations and solve (after choosing 
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proper boundary conditions) for two of the four fluctuation fields. We obtain 

1 



2a 4 (-l + 3A) (8H 2 (3ft + a 2 A) (-1 + 3A) + n 2 (-2k + a 2 A(-l + A)) ip 2 ) 
x {-16a 4 i? (3ft + a 2 A) (1 - 3\) 2 ip 



2a 4 (-l + 3A) (a 2 HQA(-l + 3A) + (-2ft + a 2 A(-l + A)) Qj ip 



D 



+ (2ft - a 2 A(-l + A)) (k 2 (ft - a 2 A) pi 2 (a 2 A?/> + 3ftV) + 2a 4 Q(l - 3A)V)] } , (4-7) 

K 2 

4a 2 (8H 2 (3ft + a 2 A) (-1 + 3A) + n 2 (-2ft + a 2 A(-l + A)) ip 2 ) 
x ^Aa 4 H(-l + 3A) (iHQ + Q) ip + 4a 4 (-1 + 3A)V^o - a 4 Qn 2 tpl 
+2H [n 2 (-ft + a 2 A) ^ 2 (a 2 A^ + 3fc>) + 2a 4 Q(-l + 3A)V'] } . 

The above solutions should be understood in momentum space where A = —ft 2 /a 2 . That is, we decompose the perturbations 
into eigenfunctions Qk{x) of the background spatial Laplacian: 



A + - )Q k (x)=0, 



(4.8) 



with eigenvalues 



ft 2 > , k = 

k 2 =£(£ + 2), k = +l 
k 2 >l, k = -l 



(4.9) 



Note that there is no singularity at the bounce point because (po ^ at the bounce time except for a measure zero set of initial 
conditions on the phase of oscillation of ipo- 



B. Quadratic Action 
Using ( 14.7b . we get a quadratic action for the two variables ip and Q: 



S 2 [ip,Q] = / dtd 3 xy^ 



where 



ip 2 + cqQ 2 + c c ipQ 

+ flpl}>i> + fQ QQ + fc iJjQ + fc Qlp + U-4,1p 2 +UJ Q Q 2 + LU C IpQ 

g = det.g y = 



1 



kr 2 \ 
4 J 



(4.10) 



(4.11) 



and the various "coefficients" (whose explicit expressions are given in Appendix IA 11 1 should be understood in momentum space. 
From Appendix IA ll we notice that 



c v Q 2 + ip 2 + c c Qip oc I V + —Q 



(4.12) 



which means that there is in fact only one dynamical degree of freedom in our system. This degree of freedom is precisely the 
Sasaki-Mukhanov l24ll combination of matter and metric perturbations, defined as 



- c = V + 



H 
— ( 



(4.13) 



which is the gauge-invariant curvature perturbation on uniform-density hypersurfaces 5 . From ( 14.131 ). we can express Q in terms 
of ifi and £, 

Q = -f(C + V0, 

Q = - I ^5 I (C + VO - Jj [C + V>J , etc. 

After plugging the above relations into ( 14. 10b . using the background equations of motion and performing many integrations by 
part, we get a new action for the two variables £ and ip: 

S 2 [C, ^} = f dtd 3 x ^ (c c C 2 + h CC + fc& + w C C 2 + u c Cip + Q^ 2 ) (4.15) 



where the various coefficients can be found in Appendix IA 21 

From (14.15b . it follows that ip is not an independent dynamical variable but rather a pure constraint, which can be solved for 
explicitly in terms of ( 

V> = wz ■ (4.16) 

Note that the coefficient Q does not vanish in our background. Hence, there is no strong coupling instability related to the 
constrained field ip. 

Plugging ( 14.161 ) into ( 14.151 ). we get an effective second-order action for a single variable £ 

S 2 [(] = J dtd 3 xV§(r( 2 + f(( + LUHLC 2 ) , (4.17) 

with 

f 2 

fr io r 

/S *"W' (4 - 18) 

The reader can verify that in the limit A = 1 and for vanishing spatial curvature the coefficient T reduces to that in the general 
relativistic theory of cosmological perturbations 
After integrating by parts, we have 

s 2 [C] ^ J dtd 3 xV$ (r c 2 - fi c 2 ) , (4. 19) 



where we made use of the definition 



fi 



(4.20) 



2 

In order to write the action in canonical form, we introduce the new variable 

u = z(, (4.21) 



5 Note that there is a singularity in the defining equation for £ at times when ipo = 0. This singularity is due to the fact that at these times the uniform density 
hypersurface becomes degenerate and hence f ceases to be a good variable to describe the fluctuations. This problem also arises during reheating in inflationary 
cosmology and in that context was studied in detail in BtZril with the conclusion that f continues through this singularity without any problem. 
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with z = a-\/2T. After changing to conformal time 77 (which is defined by dt = adrf) we have 

S 2 [u] = J di]d 3 x\/Ti^ 

where H = a' /a, h = det hij, hij is the background spatial metric without the factor a 2 , and a prime indicates the derivative 
with respect to conformal time. Note that the above result should be understood in momentum space. 
The classical equation of motion for the canonically normalized variable u is simply 

u'l+uj 2 { V ,k)u k = 0, (4.23) 

with 

u; 2 (i h k) 



As a consistency check, it can be verified that in the limit of General Relativity one obtains the usual equation of motion found 
e.g. in [22]. 



(4.22) 



a 2 n z" 




(4.24) 



V. EVOLUTION OF PERTURBATIONS DURING THE BOUNCE 



In this section we will study the evolution of fluctuations from the time of their generation early in the contracting phase until 
late times in the expanding period. We will first review the evolution of fluctuations in a matter-dominated phase of contraction 
in General Relativity. Then, we study the changes to the evolution in the matter-dominated contracting phase which arise when 
the dynamics is studied using the equations of Hofava-Lifshitz gravity. The third step is to study the dynamics of the fluctuation 
modes in the bounce phase, the transition period between matter-dominated contraction and matter-dominated expansion. Finally, 
we need to match the solutions in the bouncing phase to those in the post-bounce matter period. 



A. Einstein Gravity Analysis of the Matter-Dominated Contracting Phase 

Since we (for the sake of simplicity) model matter in terms of a scalar field, we describe the background matter in terms of a 
scalar field condensate (/Jo- 
in the limit H 2 <C m and making use of the WKB approximation, the equation of motion for ipo 

tpo + SHipo + m 2 ip = , (5.1) 

can be solved, and the solution is 

^ocm-^a-^eMijmdt), (5.2) 

and thus we see that the energy density is propotional to a~ 3 and the time average of the pressure is approximately equal to zero. 
Thus, the oscillating scalar field condensate indeed gives us a matter-dominated contracting background cosmology. 

We now turn to the description of the curvature fluctuations in the matter-dominated contracting phase, first making use of 
the perturbation equations from General Relativity. In this case, the equation of motion for the canonical fluctuation variable u 
defined in ( 14. 2U reduces to 

u'l + u; 2 ( V ,k)u k = 0, (5.3) 

with 

to 2 = k 2 - — , (5.4) 

z 
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where zaaas long as the equation of state of the background is unchanged. In a matter-dominated phase, then on scales much 
larger than the Hubble radius (where the first term on the right-hand side of ( 15.41 ) can be neglected) we have 



2 
'/ 2 



u 2 = -^. (5.5) 



If the fluctuations originate as quantum vacuum perturbations on sub-Hubble scales early in the contracting phase, then at 
Hubble radius crossing we have a vacuum power spectrum P u for u, i.e. 

u k (VH{k)) = -=, (5.6) 
V 2fc 



where r}H{k) is the conformal time when the scale k exits the Hubble radius, and hence 

Pu{m{k)) ~ k 2 . (5.7) 

To convert this vacuum spectrum into a scale-invariant one, we require a mechanism which boosts the amplitude of long wave- 
length modes relative to those of short wavelength ones. 

In an expanding universe, the amplitude of the dominant mode of u on super-Hubble lengths grows as z(rj) and hence the 
curvature fluctuation £ is constant. In contrast, in a contracting phase the dominant mode of u grows as r/^ 1 (and hence ( 
scales as z" 1 ??" 1 ). This provides exactly the boost of long wavelength modes required to turn the initial vacuum spectrum of 
fluctuations into a scale-invariant one, as can be seen as follows: 

P c (k, V ) = z- 2 ( V )P u (k lV ) = z- 2 ( ?? )(^^) 2 , (5.8) 

where in the second step we have made use of the time evolution of u. Since in the matter phase 

m {k)^k-\ (5.9) 



the factors of k in (15.8b cancel and we indeed obtain a scale-invariant spectrum. 

Another way to reach this conclusion is to consider the equation of motion of the curvature perturbation, 

C fe +3#d-^a = 0, (5.10) 

which has a solution 

ic s [l - ic s k(-q - t} B )\ ,. , . 
Cfc = A j==. — - exp[ic s k{ri - r, B )\ ■ (5.11) 

The constant factor A is determined by the initial condition. Thus the spectrum of the curvature perturbation is scale-invariant in 
the contracting phase. The initial conditions yield 



where Hb is the Hubble scale, Hb is the conformal Hubble scale, and the subscript "£?" denotes the momentum when the 
contracting phase ends. 

Since the spectrum of fluctuations is scale-invariant and the fluctuations exit the Hubble radius with an amplitude smaller than 
1, the amplitude of the fluctuations always remains perturbatively small. 



B. Hofava-Lifshitz Contracting Phase 

In this subsection we follow the fluctuation modes during the Hofava-Lifshitz contracting phase between when they exit the 
Hubble radius with a vacuum spectrum until the end of the contracting phase, when the higher derivative terms scaling as a~ 4 
in the action become important. The terms scaling as a~ 2 in the action which are induced by the spatial curvature are negligible 
throughout since we are starting the evolution in the matter-dominated contracting phase and therefore the curvature radius is 
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larger than the wavelength we considered. Thus, it is a good approximation to first solve the curvature perturbation in flat 
space (fc = 0), and in the limit A — > 1. In this limit, we can simplify the expression for lo 2 and obtain 

If A = 1, then 

u? = -% + c 2 k 2 . (5.14) 

For modes outside of the Hubble radius the k 2 term is negligible (note that Hubble radius crossing corresponds to krj = 
1). Hence, the correction terms to the mode equation in Hofava-Lifshitz gravity compared to those in the Einstein theory are 
negligible in the contracting phase. The same conclusion can be reached if we keep the leading terms due to spatial curvature. 
Assuming that the universe is closed, the action for the background becomes 



S a v = J dtd 3 x^/gN 



. q 3 , 9(1-4A) 

(l-3A)c 3 .-2 + 8(r=3iy ,31 3A 



2ttG l Aa 4 4(l-3A)a 2 8(1-3A) J 



dtd 3 x^N [ ^J' P(l)— + 0(1)- + 0(1)A] . (5.15) 



In the limit Aa 2 <C 1 which is relevant in our case when we follow modes in a phase in which the cosmological constant has 
a negligible effect, the R 2 terms dominate. (The limit Aa 2 3> 1 would correspond to a phase during which the cosmological 
constant is dominant. In this case the vacuum energy compensates the cosmological constant term in Hofava-Lifshitz gravity.) 
In our case, Eq. (l5. 1 3b can be expanded to second order of fc, yielding 



2 



UJ 



— + c 2 (3-2A)fc 2 , (5.16) 



if 2 



from which it follows that the curvature perturbation is the same as in Eq.( 15.1 11 1, except that c s = c\/3 — 2A. 

To conclude, we see that the higher derivative terms in Hofava-Lifshitz gravity lead to correction terms in the equation of 
motion for super-Hubble scale fluctuations. 



C. Hofava-Lifshitz Bouncing Phase 

The contracting phase transits into the bouncing phase when the term in the action scaling as a~ 4 becomes important. Inspec- 
tion of the Friedmann equation ([33} shows that this happens when 

T~T (^) 2 . (5.17) 
M 

During the bouncing phase the scale factor is given by ( 13.91 ) with y ~ rj^ 2 . In this case, it follows from Eq.( 14.24b that the 
mode frequency is given by 

lo 2 = [(-4 + tt)((-4 + fc 2 )K 8 /i 4 - 2kV(2ttj/(5 - 6A - 27A 2 + fc 4 (l - 4A + 3A 2 ) - 2k 2 (3 - 11A + 6A 2 )) 
+ (-5 + /fc 2 )K 4 A A i 2 )4 + (-48ttV(2 + k 2 {-l + A))(-l + 3A) 3 

+47n/ft 4 (fc 4 (l - 4A + 3A 2 ) - 9(-l + 2A + 3A 2 ) - 2k 2 {4 - 15A + 9A 2 ))A^ 2 + (-6 + /s^AV^)]/ 
[64(-3 + k 2 )y(-l + 3A)a|(K 4 (l + fc 2 (-l + A) - 3A)/i 2 + (87ry(l - 3A) 2 + k 4 (-1 - fc 2 (-l + A) + 3A)A^ 2 )a|)r/ 2 ] 

(5.18) 

where we already assumed that k = 1. Setting A = 1 yields 

2 = M + tt) (2c 2 (fc 2 - 4) + (-2c 2 (fc 2 - 6K) - 3ny) Aa 2 B ) 
^ 16 (fc 2 - 3) yrj 2 Aa 2 B ' P ' 
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and after some approximations one obtains 

(-4 + fc 2 ) (-4 + tt)kV 
256 (-3 + k 2 ) ya 2 B r] 

With this expression for the frequency, the solutions of the mode equation become 



^ 2 = • ..... ., ,, ,j • (5-20) 



( - 4+fc 2 )( _ 4+ ^ 4^2 r (_ 4 + fc 2 )( _ 4 + ^ 4„2 

C = ci(r? - rj B )^ V -^H5VT )/2 + 4( 7? _ 7?B ) (1+ V 1+ -(-3+^ )/2 . (5.21) 



If we make the reasonable "fast bounce assumption" j^rj^ 3> 1 then we find that for IR modes ( i.e. modes with k <C 1) the 
value of £ is almost unchanged between the beginning and end of the bounce phase, i.e. 

ICel^lCd, (5.22) 

where Cc(Ce) denotes the value of £ at the end of the contracting phase and beginning of the expanding phase, respectively. 

Note, in particular, that - as follows from the equation of the curvature perturbation in the bouncing phase - there is no 
singularity or instability in the solution for £. Thus, we conclude that the fluctuations pass through the bounce without singularity 
and without change in the spectrum. 



D. Expanding Phase 

In the expanding phase, the mode equation for f has two fundamental solutions. The dominant mode is constant in time on 
super-Hubble scales, the second one the decaying. Thus, the spectrum of £ on super-Hubble scales at late times is the same one 
as emerges after the bounce at the beginning of the expanding phase. 



E. Matching Condtion 

As we have seen, there are three phases in the matter bounce senario. In each phase we have derived approximate analytical 
solutions of the mode equations. All that remains is to match them correctly. The matching conditions for cosmological pertur- 
bations across a space-like slice were discussed by Hwang and Vishniac 12711 and by Deruelle and Mukhanov 12811 . These works 
show that £ and $ must be continuous across the transition surface. 

As stressed in l29ll . these matching conditions for fluctuations are only applicable if the background satisfies the continuity of 
both the induced metric and the extrinsic curvature on the matching surface . If one were to match across a singular transition 
between contraction and expansion as was done [f30h in four dimensional toy models of the Ekpyrotic l3lll scenario, then the 
background does not satisfy the matching conditions and hence the applicability of the matching conditions to the fluctuations 
is questionable. However, in a non-singular bouncing cosmology such as the one we are considering here we can apply the 
matching conditions consistently at the transition between the contracting matter phase and bounce phase, and between the 
bounce phase and the expanding matter phase. This procedure has already been applied in the case of the nonsingular mirage 
cosmology bounce of 13211 . the higher derivative gravity bounce rf33ll . and in the quintom and Lee- Wick bounces 0. 

Matching between the contacting phase and the bouncing phase implies that the spectrum of ( at the beginning of the bounce 
phase is the same as it is at the end of the contracting phase, namely scale-invariant. Since the mode functions are the same at the 
beginning and end of the bounce phase it follows that the spectrum is scale-invariant at the end of the bounce phase. Matching 
at the transition between the bounce phase and the expanding phase preserves the scale-invariance of the spectrum. Hence, we 
conclude that the spectrum of cosmological perturbations is scale-invariant at late times. 

More specifically, the values of the mode functions in the expanding phase are given by 

Ce= i^Sp exp( ^ } ' (5 - 23) 



and hence the spectrum of curvature perturbation is 



p < = = A (5 - 24) 



where c s = v3 — 2A, and Hb is the value of |U| at the end of the contracting phase (the maximal value of |if |. 
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VI. CONCLUSIONS 



We have studied the evolution of linear cosmological perturbations in a bouncing Hofava-Lifshitz cosmology. We have seen 
that at linear order in perturbation theory there are no extra dynamical degrees of freedom, the same conclusion as was reached 
in an expanding Hofava-Lifsthitz cosmology @]. 

The equations of motion for the fluctuations contain a singularity at the bounce point. We have seen, however, that the 
solutions are non-singular and thus can be smoothly extended from the contracting to the expanding phase. We have derived 
approximate solutions of the equations of motion in the contracting and bounce phases. We have seen that the extra terms in the 
Hofava-Lifshitz action have a negligible effect on the evolution of fluctuations on super-Hubble scales. Thus, an initial vacuum 
spectrum of sub-Hubble fluctuations in the far past evolves into a scale-invariant spectrum of curvature fluctuations on super- 
Hubble scales at the end of the contracting phase. Because of the smooth matching of the fluctuations between the contracting 
phase and the bounce phase, and between the bounce phase and the expanding phase, and because of the fact that modes which 
are super-Hubble at the end of the contracting phase hardly change between the beginning and end of the bounce phase, the 
scale-invariance of the spectrum of cosmological perturbations is preserved during the bounce, as initially conjectured in lfl2ll . 

We have seen that initial vacuum fluctuations lead to a power spectrum which is perturbatively small throughout the bounce 
phase. However, if one were to arrange the value of the spatial curvature and the cosmological constant such that a cyclic 
background would result, the non-trivial evolution of fluctuations on super-Hubble scales in the contracting phase would destroy 
the cyclicity of the evolution l34ll . The fluctuations would no longer be perturbatively small during the second bounce. Thus, one 
should not consider values of the parameters in the Hofava-Lifshitz action which would lead to a cyclic background. 
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APPENDIX A: VARIOUS COEFFICIENTS 
1. Coefficients in <410l 

The coefficients which appear in Eq. ( 14.101 ) are 

= 4(3fc + a 2 A) (-l + 3A)y§ 

H 8F 2 (3fc + a 2 A)(-l + 3A) + K 2 (-2fc + a 2 A(-l + A))^ ' 

= 4ff 2 (3fc + a 2 A)(-l + 3A) 

CQ 8ff 2 (3fc + a 2 A)(-l + 3A) + K 2 (-2fc + a 2 A(-l + A))<^ 2 ' 



8H (3fc + a 2 A) (-l + 3A)<p 



8H 2 (3k + a 2 A) (-1 + 3A) + n? (-2k + a 2 A(-l + A)) tp% 



/„ = - 



(A3) 



o 4 k 2 (8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) y> 2 ) 

x [(3k + a 2 A) (24a 4 H 2 (l - 3A) 2 + (3k + a 2 A) k 4 (jfe - a 2 A) ^ 2 ) + 3a 4 K 2 (-2k + a 2 A(A - 1)) (3A - , 

(A4) 

«Vo (a 2 HA(l - 3A)<ai + (2k - a 2 A(-l + A)) V) 
}Q ~ 8ff 2 (3fc + a 2 A)(-l + 3A) + fv 2 (-2fc + a 2 A(-l + A))^ 2 ' ( ) 
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fc = 



(-1 + 3A) (-a 2 A K 2 (pl + 8H (3k + a 2 A) V) 
8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) p 2 



(A6) 



2a 4 (-l + 3A) (8H 2 (3k + a 2 A) (-1 + 3A) + n 2 (-2k + a 2 A(A - 1)) tp 2 ) 

x [(3k + a 2 A) (4&a 4 H 2 (l - 3A) 2 + k 4 (2k - a 2 A(A - 1)) (fc - a 2 A) A i 2 ) + 6a 4 K 2 (-2k + a 2 A(X - 1)) (3A - l)^ 2 ,] , 

(A7) 



8a 8 (K - 3kA) 2 (8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) 
x | (3fc + a 2 A) - 1152a 8 iT 4 (l - 3A) 4 - 16a 6 H 2 K 4 (l - 3A) 2 (fc(A(-4 + 3A) - 6A) + a 2 A(A(-l + A) - A)) ^ 2 

+ (3fc + a 2 A) k 8 (2k - a 2 A(-l + A)) (fc - a 2 A) V 

- 2a 6 K 2 (-l + 3X)tpl [24a 2 i7 2 (l - 3A) 2 (-12k + a 2 A(-5 + 3A)) 

- k 4 (2k - a 2 A(-l + A)) (fc(-4A + 3AA - 6A) + a 2 A(A(-l + A) - A)) /i 2 

- 6aV (-2k + a 2 A(-l + A)) (-1 + 3A)<pg | , 

(A8) 



64H 2 (3k + a 2 A) (-1 + 3A) - 8k 2 (2k - a 2 A(-l + A)) 



9o 



- 192H 2 kA 3 g 3 - Ma 2 H 2 A 4 g 3 + 576H 2 kA 3 Xg 3 + 192a 2 H 2 A 4 Xg 3 + 12a 2 H 2 A K 2 ipl 

- 36a 2 H 2 An 2 \tpl - 16kA 3 K 2 g 3 tpl - 8a 2 A 4 k 2 g 3 tp 2 Q + 8a 2 A 4 k 2 Xg 3 pl + a 2 An 4 ip 4 
+ 8A gi (-8H 2 (3k + a 2 A) (-1 + 3A) + n 2 (2k - a 2 A(-l + A)) 

+ 8A 2 g 2 (8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) t$) 

+ 8a 2 HAn 2 Lp V' - 2Aa 2 HAK 2 Xip V + 8kn 2 {V'f + 4a 2 Ak 2 {V'f - 4a 2 An 2 X {V'f 

+ 4 (-8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (2k - a 2 A(-l + A)) tp 2 ) V"\ , 



(A9) 



2a 4 (-l + 3A) (8H 2 (3k + a 2 A) (-1 + 3A) + n 2 (-2k + a 2 A(-l + A)) cp 2 ) 
x \a 2 HA (3k + a 2 A) n 4 (-l + 3A) (fc - a 2 A) fi 2 (p 

(3k + a 2 A) (A8a 4 H 2 {\ - 3A) 2 - n 4 (2k - a 2 A(-l + A)) (fc - a 2 A) ^ 



+ 6a 4 K 2 (-2k + a 2 A(-l + A)) (-1 + 3X)(p% 



V 



(A10) 



2. Coefficients in <4.15t 



The coefficients which appear in Eq. ( 14. 15b are 



4(3fc + a 2 A) (-l + 3X)(p 2 



8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) ip 2 



(All) 
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JC a 4 F 2 (-l + 3A) (8H 2 (3k + a 2 A) (-1 + 3A) + k 2 (-2k + a 2 A(-l + A)) ip 2 ) 

x |j? (3fc + a 2 A) [24a 4 i? 2 (l - 3A) 2 + kn 4 (fc - a 2 A) pi 2 ] v? (A12) 
+ 3a 4 Hn 2 (-2fc + a 2 A(-l + A)) (-1 + 3X)tpl + 8a 4 H 2 (3k + a 2 A) (1 - 3X) 2 V 
+ a 4 K 2 (-2k + a 2 A(-l + A)) (-1 + 3A)v? 2 V } , 



A (3fc + a 2 A) k 4 (-1 + A) (-fc + a 2 A) pt 2 ip 2 



2a 2 F(-l + 3A) (SH 2 (3k + a 2 A) (-1 + 3A) + n 2 (-2k + a 2 A(-l + A)) ip 2 ) 



(A13) 



16a 8 # 3 (-l + 3A) 3 (8H 2 (fc 2 - 3fc) (-1 + 3A) + n 2 (2k + fc 2 (-l + A)) <p 2 ) 
x (-6a 8 Hn 4 (2k + fc 2 (-l + A)) (1 - 3A) 2 ^+ 

16a 4 (fc 2 - 3fc) (H - 3HX) 2 (24a 4 H 2 (l - 3A) 2 + kn 4 (fc - a 2 A) pi 2 ) (p V + 

2a 4 n 2 (-l + 3A) (8a 4 (H - 3HX) 2 (l2fc + fc 2 (-5 + 3A)) + kn 4 (2k + fc 2 (-l + A)) (fc - a 2 A) pi 2 ) tp\V' - 
4aV (2fc + fc 2 (-l + A)) (1 - 3A)VoV" + 64a 8 ff 3 (fc 2 - 3fc) (1 - 3A) 4 (V'f + 

2a 2 Hn 2 {-\ + 3X)(po (l2a 6 (H - 3H\) 2 (l2fc + fc 2 (-5 + 3A)) + 3a 2 fc K 4 (2k + fc 2 (-l + A)) (fc - a 2 A) pi 2 - 
4 (2k + fc 2 (-l + A)) (1 - 3A) 2 (-2a 4 fc 2 .gi - 2a 2 fc 4 52 + 2fc 6 .a 3 + a 6 V")) + 

Htfl ((fc 2 - 3fc) (24a 4 i7 2 (l - 3A) 2 + fc K 4 (fc - a 2 A) p 2 ) 2 + l28a 6 H 2 k 2 (fc 2 - 3k) (1 - 3A) 4 . 9 i + 
128a 4 H 2 k 4 (fc 2 - 3fc) (1 - 3A) 4 a 2 + 8a 2 (-l + 3A) 3 

(-16# 2 fc 6 (fc 2 - 3fc) (-1 + 3A).a 3 + a 6 (k 2 (2k + fc 2 (-l + A)) (V^') 2 - SH 2 (fc 2 - 3fc) (-1 + 3A)V"))) 



1 

16a 8 iF(-l + 3A) 3 (8iJ 2 (fc 2 - 3fc) (-1 + 3A) + k 2 (2k + fc 2 (-l + A)) ip 2 ) 

x (-2a 2 Hk 2 n 2 (2k + fc 2 (-l + A)) (-1 + 3A) (aV (-fc + a 2 A) pi 2 - 16(1 - 3A) 2 (a 4 9l + a 2 k 2 g 2 - fc 4 53 )) ip 4 + 
32a 4 (fc 2 - 3fc) (H - 3HX) 2 (l2a 4 H 2 {l - 3A) 2 + fc« 4 (fc - a 2 A) pi 2 - 8a 4 (l - 3A) 2 i?) <p V' + 
4a 4 K 2 (-l + 3A) (4a 4 (i? - 3HX) 2 (l8fc + fc 2 (-7 + 3A)) + 

fc K 4 (2k + fc 2 (-l + A)) (fc - a 2 A) [i 2 - 8a 4 (2k + fc 2 (-l + A)) (1 - 3A) 2 i/) tplV - 

8a 8 n 4 (2k + fc 2 (-l + A)) (1 - 3A) 2 ^y' + 128a 8 H 3 (fc 2 - 3k) (1 - 3A) 4 y' U + V') + 

Hipl (-fc 2 (fc 2 - 3fc) k 4 (fc - a 2 A) pi 2 (-16a 4 ff 2 (l - 3A) 2 + fcV(-l + A) (fc - a 2 A) pi 2 ) + 
16a 2 (-l + 3A) 3 (\m 2 k 2 (fc 2 - 3fc) (-1 + 3A) (a 4 gi + a 2 k 2 g 2 - fc 4 <? 3 ) + 

a e n 2 (2k + fc 2 (-l + A)) 4> Q V + a 6 K 2 (2k + fc 2 (-l + A)) (V^') 2 ))) 
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U * = 16a 8 tf 3 K 2 (-l + 3A) 3 

(4a 4 (l - 3A) 2 ij (24a 4 ff 3 (l - 3A) 2 + H (fc 2 - 3ft) k 4 (ft - a 2 A) ^ 2 + 2aV(-l + 3A)^> (3i?<A) - 27')) + 
(k 2 (-2a 2 #K 2 (-l + 3A) (l2a 6 (iJ - 3if A) 2 (16ft + k 2 (-7 + 5A)) + a 2 (fc 2 - 3ft) n 4 (2k + fc 2 (-l + A)) 
(-ft + a 2 A) n 2 - 8 (2ft + fc 2 (-l + A)) (fc - 3fcA) 2 (a 4 gi + a 2 k 2 g 2 - fc 4 a 3 )) ip 4 + 
6a 8 Hn 4 (2ft + fc 2 (-l + A)) (1 - 3A) 2 ^ - 4aV (2fc + fc 2 (-l + A)) (1 - 3A) 2 ^F' - 

2aV(-l + 3X)if>l (24a 4 (2ft + fc 2 (-l + A)) (H - 3H\) 2 j> a + 

(4a 4 # 2 (l - 3A) 2 (fc 2 - 6fc + 3fc 2 A) - kn 4 (2ft + fc 2 (-l + A)) (ft - a 2 A) p 2 ) V') - 

16a 4 (fc 2 - 3ft) (H - 3HX) 2 ip (24a 4 H 2 (l - 3A) 2 O + (l2a 4 £T 2 (l - 3A) 2 + kn 4 (-ft + a 2 A) /i 2 ) V'^j + 

2H 3 (fc 2 - 3ft) (-1 + 3A) (l6a 4 ijV(l - 3A) 2 (-fc 4 (-l + A) + 3fc 2 fc(-l + A) - 3fc (ft - a 2 A)) /i 2 - 

(fc 2 - 3ft) k 6 (2ft + fc 2 (-l + A)) (ft - a 2 A) 2 / + 32a 8 (-l + 3A) 3 V (<fo + V ')) + 

Hip 2 (-960a 8 H 4 (fc 2 - 3ft) (1 - 3A) 4 - 4a 4 i? 2 K 4 (l - 3A) 2 (fc 6 (-l + A) 2 - 30fc 2 (fc - a 2 A) + 

fc 4 (ft(-9 + (8 - 3A)A) + 4a 2 A) - 3fc 2 fc (ft(-9 + A) + a 2 (7 + A) A)) ^ 2 - 
(fc 2 - 3ft) k K s (ft + fc 2 (-l + A)) (fc - a 2 A)% 4 + 8a 2 (-l + 3A) 3 (l6ff 2 fc 2 (fc 2 - 3ft) (-1 + 3A) 
(aV + a 2 k 2 g 2 - fc 4 53 ) + a 6 /t 2 (2ft + fc 2 (-l + A)) 4> V + a 6 K 2 (2ft + fc 2 (-l + A)) (^') 2 )))) / 

(8H 2 (fc 2 - 3ft) (-1 + 3A) + k 2 (2ft + fc 2 (-l + A)) $)) 

These expressions simplify dramatically if one sets A = 1 and fc = 0. Setting, in addition, w = as is appropriate for the 
contracting phase w = 0, one gets 



fc = 

fc = o 



G 

36 



LOr 



■Jj r 



r] 3 n 2 aB 



6 (9 + k 2 rf 
rfn 2 a 2 B 



3fc 2 (-. 4 AA/ 2 + fffc) 
8k 2 

3fc 2 (-k 4 AM 2 + 



16k 2 



where M = fj,/a 

Inserting these coefficients into the expression for the frequency uj one obtains 



2 

n 

If the condition A = 1 is slightly relaxed and one expands to first order in A — 1 then one gets 

4V^ c2fc2 f-24 + c 2 fcV-^ 
rj 1 ) 12 \ ri A Ka 



uj 2 = -— + c 2 k 2 . (A14) 



- 2 = ( + tV 2 * 2 - 24 + ^ Tjhr ) (A - 1) • (M5) 
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APPENDIX B: USEFUL FORMULAE 
1. Preliminaries 

The connection for the perturbed metric ( 14.2b is (up to second-order): 

k 



ij'(O) 



(x J S lk + x k 5ij -x l S jk ) 



2(l + ^) 

r- j(1) = - (djipSik + dkt/jSij - diipSjk) , 

The Laplacian for the metric ( 14.2b including fluctuations is A = Am) + Am + ■ ■ • , with 



(Bl) 



A( ) = 5(0)^^- - 9(o) r «(o) a ft . 

,«j pi; 



(B2) 



where 5(o) i; , = ffij- F° r our purpose, we do not need higher-order terms such as A( 2 ) etc. Obviously, A( ) is just the background 
Laplacian (acting on scalar fields) and An) is the first-order backreaction on the Laplacian due to the fluctuation i/j. 



2. Expansions of Various Quantities 



For the metric (14.2b . we have 



£ = E tj E l] - XE 2 



3(1 — 3A) ( H 







1 - 2ip 



2(1 -3A) H 



1> 



1 - 2<ip 



AB 



ViVjBV^B -X(ABY 



(B3) 



and 



1 



1 - 2ip 



2k , . . , 38 k ^d k ^ 



(B4) 



with diipd l ip = diipdjip. The above results are exact. 

Now it is straightforward to expand various quantities up to second-order in perturbations. In this appendix we simply collect 
the final results. 



• £ = EijEv - XE 2 



Denote £ = EijE'^ — XE 2 as a shorthand. Then, to second-order in perturbations, we have 

£ (0) = 3(1 - 3A)tf 2 , 
£ (1) = -2H(1 - 3A) 3^ + A [Q) B 

£ {2) = (1 - 3A) hip 2 + 2ijjA {a) B - l2Hipi() ~ QH^A {Q) B 



(B5) 



9 2k 

(1-A) (A (0) B) 2 + ^BA (0) 5 



i? 





6fc 


#(0) 


= ^2 




4 


%) 


= ^ 


^(2) 





(B6) 



12/fc 



Ho) 1 
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IF 



(0) 



R 2 



R 2 



(i) 



(2) 



36k 2 

a 4 
48fc 
If 



a 2 A (0) V- + 3^) , 



(B7) 



/ . a2 54/c 2 , 2 27A: , A 



RijR 13 



RijR io ) 




12fc 2 


(0) 


~ a 4 


RijR io ) 




16fc 


(i) 


= 


//,,//'') 


(2) 


= 6(A ( 



(B8) 



74/c 



144fc 2 
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